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Coupling between whistler waves and slow-mode solitary waves 
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The interplay between electron-scale and ion-scale phenomena is of general interest for both 
laboratory and space plasma physics. In this paper we investigate the linear coupling between 
whistler waves and slow magnetosonic solitons through two- fluid numerical simulations. Whistler 
waves can be trapped in the presence of inhomogeneous external fields such as a density hump or 
hole where they can propagate for times much longer than their characteristic time scale, as shown 
by laboratory experiments and space measurements. Space measurements have detected whistler 
waves also in correspondence to magnetic holes, i.e., to density humps with magnetic field minima 
extending on ion-scales. This raises the interesting question of how ion-scale structures can couple 
to whistler waves. Slow magnetosonic solitons share some of the main features of a magnetic hole. 
Using the ducting properties of an inhomogeneous plasma as a guide, we present a numerical study 
of whistler waves that are trapped and transported inside propagating slow magnetosonic solitons. 
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I. INTRODUCTION AND AIMS 

The study of the interplay between phenomena at ion and at electron dynamical scales is crucial in order to 
understand the physical mechanisms of basic plasma processes such as energy dissipation and particle acceleration 
and heating in space and laboratory plasmas. Presently multi-point satellite measurements such as the Cluster 
satellites, by combining space and time measurements, make it possible to investigate simultaneously electron-scale 
and ion-scale phenomena and to inspect stationary and propagating magnetic and density structures. 

An important example of coupling between ion-scales and electron-scales is provided by the correlation that has 
been found in space between whistler waves, which occur on the electron-scales, and magnetic field depressions 
associated to density humps that have typical scale lengths of the order of the ion-scales, usually interpreted as 
non-propagating mirror mode structures [IHSl. On the other hand, whistler waves can also interact with slowly 
propagating Magnetohydrodynamic structures involving both density and magnetic field modulations. If the whistler 
waves become trapped inside such structures, the problem arises how low frequency nonlinear modes can act as carriers 
for higher frequency waves. 

Whistler waves are electromagnetic right-handed polarized waves that propagate nearly parallel to the ambient 
magnetic field at frequencies in the range between the ion and the electron cyclotron frequency [6] . A known property 
of whistlers is that in the presence of plasma inhomogeneities, such as magnetic field aligned tubes of density en- 
hancements or depletions, their energy can be guided for long times without being dispersed [THTQ]. Examples of such 
ducted propagation has been found in satellite observations in the Earth's magnetosphere pTtfTS] and in laboratory 
plasmas [14 . 

Magnetohydrodynamic waves also have been studied for decades both in space and laboratory plasmas. Recently, 
growing attention has been devoted to the study of oblique slow magnetosonic solitary waves, namely of coherent 
structures at the ion-scale characterized by a density hump and a magnetic field depression p!5Hl8] . 
Despite many studies of both magnetosonic and whistler waves, their coupling has not been investigated in detail 
and is yet poorly understood. In this paper we investigate the interaction between slow magnetosonic solitons and 
whistler waves by using a two-fiuid model. Using the ducting properties of an inhomogeneous plasma as a guide, 
we propose a new mechanism of ducting and transport of whistler waves arising from a linear coupling with slow 
type magnetosonic solitons. Such a mechanism could explain spacecraft observations in the Earth's magnetosphere, 
namely the recurrent detection of whistler waves correlated to magnetic holes. 

This paper is organized as follows: in Section II we present the two main problems of interest separately: the 
whistler waves trapping and the slow mode solitary wave. In Section II A we treat the trapping of whistlers by an 
inhomogeneous magnetized plasma. The aim is to shed light on the ducting process using an equilibrium configuration 
simpler than the soliton solution. In this way we can obtain quantitative conditions to be used to estimate the ducting 
conditions when the inhomogeneity is provided by the magnetosonic soliton. In Section II B we introduce the slow 
mode solitary wave solutions of the two-fiuid equations and discuss the role that these configurations can have 
in trapping whistler waves, in Section III we describe the system of equations of the two-fiuid model used in our 
simulations, the initial conditions and the simulation parameters. In Section IV we investigate the trapping of whistler 
waves by slow mode solitons numerically and present the results. Finally, conclusions are discussed in Section V. 

II. THEORETICAL BACKGROUND 

A. Whistler Avave trapping in magnetic and density ducts 

It has been shown jTHTO] that whistler waves propagating in a cold magnetized plasma in the presence of inhomo- 
geneities axially symmetric and transverse to the ambient magnetic field of the plasma density can be channelled by 
these inhomogeneities {density duct). The variation of the index of refraction caused by the density inhomogeneity 
makes the wave trajectory bend such that the average propagation is along the duct direction. As a consequence, the 
wave becomes trapped by the duct. 

In this Section we extend these results to the case where both the equilibrium plasma density and magnetic field 
are inhomogeneous in the plane perpendicular to the direction of the magnetic field. 

The aim of this analysis is to study the problem of whistler wave trapping with a "simple" model. It will be used 
later as a reference for a more complex one, adopted in the numerical simulations, where the self consistent fields of 
the soliton are viewed by the whistler as perturbations of the ambient equilibrium. In order to mimic the configuration 
of interest (the slow mode soliton), we consider a plasma equilibrium characterized by a density hump and a magnetic 
field minimum, the magnetic hole. 

Model - We want to study the trapping properties of the equilibrium magnetic field inhomogeneities in the framework 
of a fiuid model. Therefore we consider the regime Vth,e <^ "^ph {vth,e and Vph being the electron thermal velocity 



and the whistler phase velocity respectively), and we use the cold dielectric tensor of a magnetized plasma, e. This 
simplification is convenient because the cold dielectric tensor e includes the basic effects of trapping thanks to its 
dependence on the density and the magnetic field strength (through the electron cyclotron frequency uOce)- Assuming, 
for the sake of illustration, a two dimensional spatial configuration and taking the gradients along the magnetic field 
lines to be negligible on the scale of the whistler wave length, the plasma can be represented in a slab geometry, 
with density and magnetic field gradients perpendicular to the magnetic field direction. For the sake of clarity, we 
define "magnetic hole" a fluid equilibrium with a density hump and a magnetic field depression perpendicular to the 
equilibrium magnetic field, as distinguished from a "density duct" that only has a perpendicular density inhomogeneity. 
This model allows us to highlight the basic mechanism of the whistler wave trapping and to obtain a quantitative 
estimate of the parameters to be used in the simulations, such as the value of the whistler frequency and the angle of 
propagation. We will show that a magnetic hole requires less strict conditions on the whistler wave parameters than 
a density duct in order to trap whistler waves. 

Let us consider a bump-like density profile in the x direction, perpendicular to the background magnetic field 
^y{x) which is taken to be directed along y. The background magnetic field has a minimum in correspondence to 
the density hump. A sketch of this configuration is illustrated in Fig. IT] the wave propagates in the (x, y) plane and 
is localized inside the magnetic hole. As shown in Ref. j8j, the set of the two- fluid equations for a magnetized plasma 
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Figure 1: Schematic representation of the trapping of whistler waves: the wave propagates in the plane [x^y) containing the 

magnetic field and the inhomogeneity direction and is localized inside the density hump. In the case of a warm plasma, a 

magnetic field depression of the form By[x) is associated to the density enhancement. 

can be arranged so as to obtain the following system of equations for the electric and magnetic fields, where a time 
dependence of the form exp(— zcjt) has been assumed: 



V(V-E)- V^E: 



-(e-E) 



(1) 



V • (£ • E) = 



(2) 



-i-V X E. 

UJ 



(3) 



Looking for spatial solutions of the form A{x) e-xjp{ik\\y) ^ where k\\ = ky is the parallel wave vector, equation (|T|) can 



be reduced to two coupled differential equations of second order for the electric field components Ex and Ez^ while the 
third component Ey is obtained from equation Q. By rescaling the variable w = x/L, where L is the typical length 
of the large scale inhomogeneity such that {c/uj)/L <C 1, the two coupled equations for Ex and Ez can be solved by 
means of the WKB approximation [8] . We impose solutions of the form 



A{w) = B{w) exp I iLuj/c I 



q{w')dw' 



and retain only the terms to lowest order in {c/uj)/L. Note that such a choice excludes the ducting mechanism 
described in [19 where the whistler energy propagates in the form of "surface waves" with frequency below the 
characteristic frequency of bulk whistler waves. The perpendicular wave vector k^ = {uj/c)q = kx, for a given 
parallel wave vector k\\ and frequency cj, must satisfy the whistler dispersion relation of a homogeneous plasma (in 



ion normalized units, see also Section III) obtained in the limit ujci < uj < ujce <^ uOp 



kl±{x) 



2dl {uj/ujce{x)) 



k\\d^ 



UJ 



OJce{x)^ 



) - ^"(^) (^) ± '^'^^ f ^^i - ^"(^) (^) • (4) 



Note that, using ion normalized units, d^ = {me/rrip) and Uce = Brrip/me. 

The whistler dispersion relation as expressed by equation Q for the complex variable k± shows that, for fixed 
/c|| and cj, there are two "branches" of the perpendicular wave vector corresponding to the plus and minus sign, 
respectively (the "upper" and the "lower" branch). For a solution given by a real k±^ we get a propagating whistler 
wave while for an imaginary k±^ we get a purely evanescent (non propagating) whistler wave. The transition within a 
given branch from real to imaginary values is at the basis of the wave trapping, and is determined by the local values 
of the density n and of the parameter uj/ujce- 

The general WKB solution is given by a linear combination of the four wave solutions corresponding to the four 
possible wave vectors ±kj_^±. Near the critical points, where k']_ ^ = or /c^ _ = k'j_ _^^ the WKB approximation 
ceases to be valid. An analytical continuation of the solution in the complex x plane around these points is therefore 
necessary in order to extend the solution to all its domain of validity. The complex x plane is divided into different 
portions by the so called Stokes and anti-Stokes lines which radiate out from the critical points [20]. When a Stokes 
line radiating from a turning point of a given branch is crossed, the two solutions of the given branch, say, the ones 
corresponding -\-k_\_^- and to —k^-^ are coupled, corresponding to the reflection of a given wave when approaching 
a turning point. In an similar way, when a Stokes line radiating from a conversion point is crossed, the coupling 
between the two branches, for example +/cx,+ and +/cx,-, occurs. The coupling between the two branches leads, for 
instance, to the leakage of a propagating wave, again approaching a turning point. Because of the coupling with the 
other branch at the conversion point, a fraction of the wave energy continues to propagate past the turning point [8]. 
The coefficients of reflection or conversion are exponentially small if the critical points are far from the real axis. 

In this Section we consider only the trapping of whistler modes due to the presence of turning points on the real 
axis, where k\ j^ =0. For simplicity the conversion between different branches will be neglected here, assuming that 
the square root in equation Q does not vanish on the real axis. Moreover, in order to look for solutions given by 
propagating waves or evanescent waves, we assume that k\ ^ is real, which means that the perpendicular wave vector 
is real (i.e. /c^ ^ > 0) or imaginary (i.e. k\ j^ < 0). These assumptions correspond to the following condition: 

^>'^{^S- (5) 

From equations (|4|-([5| it follows that the upper branch cannot be trapped in a magnetic hole (neither in a density 
duct). Indeed, k'^ ^ is everywhere positive for frequencies uj/uce < 1/2, and thus the wave propagates in all regions. 
If cj/cjce > 1/2 then k'^^_^ is positive when /cndg > n/[ijOce/^ — 1]- If this condition is satisfied inside the magnetic hole, 
it is satisfied outside the magnetic hole as well, since n{x) has lower values outside than inside the magnetic hole, and 
vice versa for the function (Jce{x). Then also in this range of frequencies the wave propagates in all regions. We can 
therefore focus only on the lower branch k±^-. 

With the same reasoning as above, we see that the lower branch can be trapped in a magnetic hole (or in a 
density duct) only for frequencies uj/uce < 1/2. Indeed, for frequencies u/uce > 1/2 the perpendicular wave vector 
corresponding to the lower branch is imaginary everywhere while for oj/oJce < 1/2, the perpendicular wave vector is 
real when 

^I < ^7^. (6) 

II UJce/^ - 1 



while it is imaginary when 



> .J^, (7) 

II UJce/(^ - 1 



To summarize, for frequencies u/uce < 1/2, trapping is possible in a magnetic hole (and in a density duct) if the 
parallel wave vector satisfies equation (Isl) everywhere (which means that k± is either imaginary or real), and safisfies 
equation ([6| inside the magnetic hole, giving a propagating wave, and equation ^ outside the magnetic hole, giving 
an evanescent wave. Since we are interested in whistler modes trapped into the magnetic hole, from now on we 
consider only the lower branch in the frequency range u/uce < 1/2 and drop the subscript "— ". For the sake of 
clarity we define riin and Uout the density calculated at the center of the magnetic hole (thus in correspondence to 
the minimum of the magnetic field and to the density maximum) and outside the magnetic hole (where the medium 
is homogeneous), respectively. By analogy we define the frequencies ujI^ and cj^^^, the electron cyclotron frequency 
calculated at the center and outside the magnetic hole, respectively. In this way, the trapping condition for the lower 
branch in a magnetic hole can be written as follows: 



where 



and 



hnf{^) 



max 









1/2 



^out 



Vw-1_ 



1/2 



M 



di w^/w-r 



(9) 



(10) 



In Fig. pi we show a graphical representation of the portions in the parameter space (cj/cj^g, k\\) corresponding to 
real values of k^ , calculated at the center (solid lines) and outside (dashed lines) the channel provided by the magnetic 
hole or the density duct. Red and black lines correspond to the right-hand-side of equation ([5| and equation (l6|, 
respectively. The left panel corresponds to a plasma equilibrium with a magnetic hole {AB/B =\B'^'^ — ^ontj^^onf^ 
0.3 and An/n = \nin — nout\/^out = 0.37) and the right panel to a density duct with the same density inhomogeneity 
than the magnetic hole {AB/B = and An/n = 0.37). Referring to Fig. pi left panel, the points {uj/ujI^^ ku) lying in 



the portion A+B and B+c correspond to a propagating wave in the region mside and outside the channel, respectively. 
The intersection B of these two regions corresponds to the untrapped modes, as they propagate both inside and outside 
the channel. The trapped modes are those corresponding to the portion A, where k± is real inside and imaginary 
outside the channel. The maximum angle Omax{^) of trapped modes for a given frequency is determined by kinf^ and 
by the corresponding k±{uj, kinf)'. 



Omax{oj) = arctan 



k±{uj,kinf) 



fvi 



inf 



(11) 



A comparison between the magnetic hole, left panel in Fig. [2J and the density duct, right panel in Fig. [2J shows 
that the presence of magnetic variations (magnetic hole) leads to less strict trapping conditions. Indeed, for an equal 
density variation, the portion of trapped modes in a channel provided by both density and magnetic inhomogeneities 
is larger than in a channel formed only by a density inhomogeneity. In addition, the maximum angle of trapping (not 
shown here) results to be higher. 





Figure 2: Plot of the curves which bound the region in the parameter space (cj/cjce 5 ^||) corresponding to real values of k±^-. 
Solid and dashed lines are calculated in correspondence to the region inside and outside the channel, respectively. Red (straight) 
and black (curved) lines correspond to the right-hand-side of equation ([5]) and equation (pi, respectively. The portion A+B and 
B+c correspond to real values of k± inside and outside the inhomogeneous region, respectively. The portion A corresponds to 
the trapped modes. Left panel: plasma equilibrium with a magnetic hole; right panel: plasma equilibrium with a density hump 

on a uniform magnetic field. 



B. Slow magnetosonic solitons acting as carriers of w^histler w^ave energy 



We are interested in slow mode solitary waves because, as it will be explained in the following, they can trap whistler 
waves, similarly to magnetic holes. In addition, since solitary waves propagate almost unchanged, they provide an 
efficient channel that not only confines but also transports whistler energy at the typical soliton propagation speed. 



Magnetosonic solitons are nonlinear one dimensional perturbations propagating in a warm plasma, obliquely to 
the equilibrium magnetic field [l'5l HSl Hi] . Solitons are characterized by magnetic field and density perturbations in 
phase (fast solitons) or in opposition of phase (slow solitons). Solitary waves propagate with a constant profile and 
arise when the non linear terms are balanced by the dispersion terms. In a two-fluid model the required dispersion 
which gives rise to magnetosonic solitons is given by the Hall term and the electron inertia. Nonetheless, for non 
perpendicular propagations, the Hall term dominates the dispersion and the typical scales of solitons are ^ di. It can 
be shown that at some level of approximation the system of two-fluid equations can be reduced to a Kortveg de Vries 
equation [21 , which has solitary wave solutions. 

Let us consider a soliton moving in the positive x direction in a homogeneous magnetized plasma at rest, with 
equilibrium quantities defined as follows: 



B = Bo = {Bqx, Boy^ 0) Ui 



(0, 0, 0) 



no 



Pi. 



Po, 



where u^, e and Pi^ e are the ion and electron velocity and pressure, respectively, n the density and B the magnetic 
field. The angle of propagation of the soliton is defined as the angle between the direction of propagation, namely the 
X direction, and the equilibrium magnetic field Bq. For future convenience, in order to define the soliton direction of 
propagation, we will use the angle (/?o, taken as the complementary angle to the angle of propagation, thus defined as 
(fo = diTctdin{Box / Boy) . An explicit solitary solution of the two-fluid system of equations is obtained in Ref. [21] in 
the weakly non linear approximation. Both fast and slow solitons are found. Here we consider only the slow mode, 
as this mode has density and magnetic field in opposition of phase. The slow mode solitary solution depends on 
the amplitude of the perturbation A, on the angle (po^ on the temperature and on the Alfven speed which, in the 
homogeneous region, is equal to one in our units. The effects of the temperature enter through the sound speed Cg 
and the slow magnetosonic phase speed Vpo^ and determines the velocity of propagation of the soliton and its width. 
Below, equations (12)-(13) represent the total magnetic field and the total density of the plasma, defined as the values 
of the homogeneous equilibrium plus the fluctuations associated to the soliton, B^ot = Bq + Bgo/, ritot = no -\- Usoi- 



ritot = 1 -^risoh 



Bx 



sirnpo, 



B 



'y,tot = COSipo + 



{vf 



pO 



I) 



COS(/?o 



'^soli 



(12) 



where 



rtsoi oc 



cosh 



12fi{ifo,Po) ' 



(13) 



The other plasma quantities are given explicitly in Appendix [A] including the expression for the function jii{(fo^Po). 
Quantities are normalized to asymptotic equilibrium values outside the soliton. 

According to this theoretical analysis, the propagation speed of the soliton is Vq = Vpo + A/3 and the typical width 
is ^ ~ 2^yl2jii/A. The function /i, which determines the width of the soliton, is a growing function of the temperature, 
ranging from values smaller than, or of the order of, di to values much greater than di. Note that for the slow mode 
'^po < Cg. Thus from the last of equations ([l2|, which defines By^^^ it follows that the magnetic field perturbation is 
in opposition of phase with the density permrbation. The analytical solution for the slow soliton is valid as long as 
the propagation is not parallel (ipo = 7r/2) in which case ja equals zero (if Cg < 1) or infinity (if Cs > 1) [21 . 

To summarize, the main features of slow mode solitons are: they carry a density hump perturbation associated 
to a magnetic field depletion and propagate obliquely with respect to the background equilibrium magnetic field at 
speeds which are much smaller than that of whistler waves (greater than unity). This inhomogeneous system, which 
can be represented by an oblique solitary perturbation moving in a homogeneous plasma at rest, is more complicated 
than the magnetic hole discussed previously, which instead has purely perpendicular gradients with respect to the 
magnetic field. Nevertheless, as a first approximation, it is possible to consider the soliton perturbation superposed 
to the background equilibrium as a local and instantaneous magnetic hole for whistlers that are injected inside the 
soliton. 



III. MODEL EQUATIONS, INITIAL CONDITIONS AND PARAMETERS 



In order to describe numerically the trapping and transport of whistlers by solitary waves in a magnetized plasma, 
we use a quasi neutral adiabatic two-fluid model. The set of equations for ions and electrons, with labels i and e 



respectively, are normalized using as characteristic quantities the ion mass density nrrii, the Alfven velocity Va = 
B/y/ATinmi and the collisionless ion skin depth di = c/upi. With this choice d^ = rrie/mi and ujce = B/d^. The 
two-fluid model equations, in dimensionless form, are: 

dn/dt + V • (nU) = (14) 

dSe,^/dt + V • (5e,.Ue,.) = 0, 5e,i = Pe,^n^-^ (15) 

d{nV)/dt = -V • [n{uiUi + d^UgUe) + (Pe + Pi + ^V^ - BB)] . (16) 

(1 - d^V^)E = -Ue X B - (1/nWPe- 

d^{ui X B - (l/n)VPi + (l/n)V • [n(uiUi - UeUe)]}, 

Ue = U - j/n , Ui = U + 6^eJ/^ 5 U = U^ + dgUe, (18) 

V X E = -dB/dt , j = V X B + zf/K (19) 

In these equations E and B are the electric and magnetic field, u^^g is the ion (or electron) velocity, j the current, 
jl^^z is an external forcing current used to inject whistlers, P^^g the pressure of ions and electrons F is the adiabatic 
index. 

The dimensions of the simulation box, L^ and L^, and the resolution of the grid, dx and dy^ are chosen in order to 
find a compromise between the different time and length scales at play. Ly is chosen in order to let the whistler wave 
train propagate for several tens of cj"^, without reaching the boundaries. L^ is chosen in order to contain the soliton 
which is wider or of the order of di and moves at a speed V ~ 0.1. Finally, the mesh size must resolve the whistler 
wavelength. In Table [T| we report the parameters of the simulation box. The ion to electron mass ratio is fixed to 
rrip/me = 100. 
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Table I: Parameters of the injected whistler mode, of the simulation box and the characteristic width of the soliton £. 

As initial condition, we consider a slow mode solitary wave centered in the simulation domain and superposed to 
a homogeneous magnetized plasma at rest. Oblique whistlers, as explained in the following, are injected artificially 
in the simulation box, during the initial phase, in correspondence to the soliton. In order to do this we make use of 
an oscillating forcing current lasting over a characteristic time r. In Fig. [3] we show a schematic view of the system. 
The dashed lines indicate the region filled by the soliton moving in the positive x direction with velocity V; BJ^^ is 
the total magnetic field at the center of the soliton, forming an angle (p with the y axis, and k is the whistler wave 
vector. In particular, when whistlers are generated inside the soliton, the subscripts "±"and "||" of the wave vector 
refer to the total magnetic field BJ^^ at the center of the soliton. Outside the soliton the total magnetic field reduces 
to the equilibrium magnetic field Bq forming an angle cpo with the y axis. 

The analytical form of the fields representing the soliton superposed to the homogeneous equilibrium at the initial 
time are given in Appendix |A] The parameters of the soliton, listed in Table [ill ^^^ chosen in order to have a narrow 
soliton with a width i of the order of the ion skin depth (simulations 1 to 6) or a wider soliton with i of the order of 
several ion skin depths (simulations 7 to 9). 




V 




Figure 3: Schematic view of the system. The dashed hnes indicate the region filled by the soliton moving in the positive x 
direction with velocity V] BJ^t is the total magnetic field at the center of the soliton, forming an angle Lp with the y axis, and k 
is the whistler wave vector, which forms an angle with the magnetic field BJ^^. Outside the soliton the equilibrium magnetic 

field Bo forms an angle lpq with the y axis. 



To our knowledge, the numerical stability of the approximate solutions discussed in Section |IIB| has never been 
investigated. Therefore, before considering the full problem including the injection of the whistlers, we tested the 
stability of the soliton solutions numerically. Our results show that they are well stable in the range of propagation 
angles ip^ = 0.57 — 0.17 and of typical variations with respect to the equilibrium Usoi ^ 0.02 — 0.8, \By^soi\/Boy ^ 
0.01 — 0.5. They propagate at the expected velocity maintaining almost unchanged their initial profile over times 
t ~ 1000, until they exit from the simulation box. 

In Fig. [4] we show two examples of quasi perpendicular {(po = 0.17) magnetosonic solitons propagating along the 
X axis at three different times up to t ~ 1000. The red lines represent the density profile Utot and the black lines 
the magnetic field By^tot- In the left panel we represent a narrow, strong amplitude soliton {An/n ^ 0.9) and in the 
right panel a wider and weaker soliton {An/n ~ 0.2). Notice that the initial soliton profile slightly modify during 
the temporal evolution, especially for large amplitude solitons, since the analytical profile is not an exact solution of 
the two-fluid system. For the sake of clarity, we indicate with a subscript "tot" the quantities resulting from the sum 
of the homogeneous background equilibrium plus the soliton perturbations. These large scale variation fields can be 
considered as the "inhomogeneous equilibrium" supporting the whistler waves. 
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Figure 4: Examples of slow magnetosonic solitons at three different times: the black lines (depletions) represent the magnetic 
field By ^ tot and the red lines (humps) the density ntot- In the left panel we represent a narrow, strong amplitude soliton 
(An/n ^ 0.9) and in the right panel a wider and weaker soliton (An/n ^ 0.2). Notice that the initial soliton profile slightly 
modify during the temporal evolution, especially for large amplitude solitons, since the analytical profile implemented as initial 

condition is not an exact solution of the two- fluid system. 



Finally, we generate small amplitude oblique whistlers with frequency ujq and propagating at an angle with respect 



Sim. 1-6 0.073 0.05 0.17 2 0.85 1.87 0.2 0.1 86 
Sim. 7-9 0.034 0.5 0.17 13 0.7 1.18 0.24 0.1 72 

Table II: Parameters of the solitons. 

to the total magnetic field during a characteristic time scale r by an external forcing corresponding to a current along 
the z axis, J|^^ (see equation (19)). The parameters of the forcing and thus of the injected whistlers, listed in Table[l[ 



are chosen making use of the trapping conditions discussed in Section II A In particular, we use as reference the 



maximum trapping angle defined by equation (11). The maximum angle of trapping at a given frequency cjq, that we 
define as 0":^^^^ allows us to choose the frequency cjq and the angle properly in order to inject a specified whistler 
mode that we expect to be trapped or not. 

IV. TRAPPING OF WHISTLER WAVES BY SLOW MAGNETOSONIC SOLITONS: NUMERICAL 

RESULTS 

In this section we show, by means of numerical simulations, that whistlers can be trapped and transported away 
by a slow magnetosonic soliton. Even if a slow soliton propagating in a homogeneous magnetized plasma is more 
complicated than the so called magnetic hole, as a first approximation the same properties of whistler ducting apply. 



and the trapping conditions found for the magnetic hole discussed in Section [nA] are therefore a good reference when 
asking which whistler modes can be trapped by the soliton. 

We have investigated the slow magnetosonic ducted and unducted regime of whistler modes by varying the typical 
width of the soliton. Here we report two different typical cases: a narrow soliton of width £ ^ 2 < k~^ and a wider 
soliton of width ^ ~ 13 > k~^ (in units of d^), where k is the whistler wave vector estimated for a given frequency 



and propagation angle from the two-fiuid cold dispersion relation (see equation (Bl) in Appendix |B|). A list of the 
parameters used in the simulations for the "small" and "large" soliton are listed in Table [TTj The wider soliton has 
a weaker density hump but a stronger magnetic field depression than the narrow one. The injected whistler modes 
fiuctuate at low frequencies {ujq <C 0.1cj*g) or high frequencies {ujq ~ O.lcj^g) with different angles of propagation 
ranging from <^ ^max to ^ > O';^^^. In the following, we focus on two simulations, namely Sim. 1 for the narrow 
soliton and Sim. 7 for the wide soliton, to show the trapping of whistlers. 

In these simulations the injected whistlers have frequency ujq ~ 0.03 cjjg and ujq ~ 0.04 cj^^ respectively. They are 
injected along the y axis, slightly oblique with respect to the local total magnetic field, forming an angle = —0.198 
and = —0.24, respectively, then satisfying \6\ <C 0':^^^. The forcing current oscillates at the center of the simulation 
domain and switches off exponentially on a characteristic time shorter with respect to that of the soliton propagation. 
In this way, the forcing generates two finite size wave packets in the (x, y) plane propagating away from the source 
region in the two opposite directions, namely in the positive y direction (upward) and in the negative y direction 
(downward). The two wave packets propagate upward or downward, respectively, and remain spatially confined along 
the inhomogeneous x direction in correspondence to the soliton, following its displacement along x. This is shown 
in Fig. |5] by the contour plots of the x component of the magnetic field bx of the whistler waves in the simulation 
domain when the current has switched off and the wave packets are well developed. The profile of the soliton is 
represented (not in scale) by black lines and the dashed line corresponds to the soliton at time t = 0. The left panel 
represents the two wave packets at time t — 100, for the narrow soliton. In the middle panel we show the propagation 
of the same wave packets as injected in Sim. 1, but in a homogeneous equilibrium, i.e. with B^y = By^^^^ = 0.85 and 
no = nj^^ = 1.87. We see that in the absence of the soliton the injected wave packets spread out during propagation. 
Finally, the right panel shows the wave packets at time t = 60, for the wide soliton. To summarize, our simulations 
provide evidence that the waves, initially injected inside the soliton structures, propagate along the duct provided by 
the soliton, upward or downward, advected at the same time in the perpendicular x direction by the soliton. The 
whistlers are thus confined and transported by the slow soliton over times much larger than their typical time scale. 

We consider now two simulations, Sim. 3 and Sim. 6 where we inject highly oblique whistlers with > 0^^^. In this 
case the waves escape outside the solitons. An example is shown in Fig. |6J left panel, where we show the contour of 
hx and the profile of the soliton at time t = 16 in the case of Sim. 6. Here the injected whistlers are at high frequency 
(cjo ~ O.lcj^g) and at an angle of propagation = 1.3. 

These results are in good agreement with the ducting theory. However, the model used in our numerical study is far 
richer than the reference model of the magnetic hole (Section |II A ), and there are important effects that can modify 
the trapping conditions. 

First of all, even if we neglect the displacement of the soliton, there is a finite perturbation in the plasma velocity of 
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Figure 5: Contour plot of the x component of the whistler magnetic field bx- The profile of the soliton is represented by the 
black lines. The dashed line represents the soliton at time t — 0. Left panel: trapping of whistlers in Sim. 1 at time t — 100. 
Middle panel: the whistler wave packets in the entire space domain at t = 30 as in Sim. 1 but in the absence of the soliton. 

Right panel: trapping of whistlers in Sim. 7 at time t = 60. 



the form \J{x) (see equations ( |A4[ )-(A5) in Appendix [A|). The presence of the fluid velocity introduces an asymmetry 
in the system, due to the Doppler effect, between wave packets propagating upward and downward, ky > and ky < 0, 
respectively. Second, there are gradients along the total magnetic field, which can drive whistlers outside the soliton 
even if trapping conditions are satisfied. In the case of the magnetic hole /cy and uj are fixed quantities, while for the 
soliton UJ and ky are constant but the parallel wave vector varies as the whistler propagates inside the soliton. As a 
consequen ce, w hile the whistler propagates towards the edge of the soliton, k\\ can approach the value kinf (as defined 
in Section II A, see also Fig. [2I left panel) thus allowing the whistler to become untrapped. An example is given in 
Fig. |6| right panel, which refers to Sim. 9. In this simulation only the lower wave packet is trapped while the upper 
wave packet is guided outside the soliton. An interpretation of Sim. 9 can be given in terms of geometrical optics. 
Since the soliton moves along the x axis at a speed V ~ 0.1 much smaller than the whistler phase velocity (greater 
than unity), as a first approximation we neglect the displacement of the soliton. Because of the Doppler shift, the 
frequency ujq measured in the simulation is given by 



cJo(k, x) = a;(k, x) + k • U(x), 



(20) 



wher e a;( k, x) is the whistler two-fluid dispersion relation in a plasma at rest obtained in the cold limit (see equa- 
tion (Bl ) in Appendix [B|). The dispersion relation a;(k, x) is given in terms of ky and k^ and the density and magnetic 
field profiles are given by n = 1 + Ugoi and B = Bq + ^soh respectively. In the framework of the geometrical optics, 
the contours of uq in the plane (k^^x) for fixed ky represent the orbits of the whistler wave packet. The solution of 
the Hamiltonian system 

duj • , . duj ... 

gives the evolution of the wave vector and the trajectory of the wave packet. In Fig. [7| we show the contours of ujq 
as defined in equation (20) obtained using the soliton profile of Sim. 9. The wave vector ky can be estimated from 
the forcing frequency and injection angle taking into account the doppler shift, giving ky ^ 1.3 and ky ~ —1.7. The 
contours in Fig. [7| show that the orbit corresponding to ct;o = 3 is open for the wave packet propagating upward while 
it is closed for the wave packet propagating downward. A Fourier analysis of the x component of the magnetic field 
in Sim. 9 confirms that the wave vectors with \ky\ ~ 1.7 are trapped inside the soliton. Similar results are obtained 
for Sim. 2 and Sim. 8. 

However, because of the movement along the x axis, the soliton behaves as a "moving mirror" thus causing the 
frequency of the injected whistler to change with time. We qualitatively estimate the shifted frequency after the first 
reflection at the soliton edge uJq = cjq — 2k • V. The change in frequency could cause the wave to become evanescent. 
An example is given by Sims. 4 and 5 in a high frequency whistler regime {uj ~ 0.1cj*g) and using a narrow soliton. 
In these simulations only the upper wave packet is trapped, as expected, while after the first reflection at the left 
boundary of the soliton, the lower wave packet becomes evanescent. Fig. |6) middle panel, refers to Sim. 4 and shows 
the contour of bx and the profiles of the soliton at time t = 30 (solid line) and t = (dashed line). In this simulation 



the lower wave packet has a wave vector ky estimated to ky 



-5 (in agreement with the Fourier spectrum of the 
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simulation results). The solution of the Hamiltonian equations for the wave packet with initial conditions x{0) = 0, 
kx{0) = gives wave packet reflected at nearly dx ~ —1.5 from the center of the soliton with kx ~ 25 (in agreement 
with the small scales which form in the x direction at the edge of the soliton). In this point the reflected frequency is 
estimated as coi^ ^ ujo — 2kxV ~ 3.5, which is below the minimum frequency calculated in correspondence to the edge 
5, explaining why the lower wave packet does not propagate after reflection. 
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Figure 6: Contour plot of the x component of the whistler magnetic field bx^ with the profile of the soliton represented by the 
black lines. The dashed line corresponds to the soliton at t = 0. Left panel: Sim. 6 at time t = 16, corresponding to a high 
frequency, highly oblique whistler that escapes from the narrow soliton. Middle panel: Sim. 4 at time t = 30, corresponding 
to a high frequency whistler that is trapped only in the upward direction while the downward wave packet becomes evanescent 
after one reflection. Right panel: Sim. 9 at time t = 60, corresponding to an upward whistler that escapes outside the soliton 
while the downward wave packet is trapped. The periodic boundary conditions cause the waves approaching the upper (lower) 
boundary of the simulation box to appear in the lower (upper) boundary. 





Figure 7: Contours of the whistler dispersion relation in the rest frame of the soliton, ujq{\^^x) — C(;(k,x) + k • U(x), in 
the (kx,x) plane with the parameters of Sim. 9. The left panel refers to the upward wave packet, which has ky ^ 1.3, and the 
right panel to the downward wave packet, which has ky ^ —1.7. Whistler wave packets with a given ky evolve moving along 
the orbits at constant frequency. In Sim. 9 the frequency is cjq = 3 that corresponds to an open orbit for the upper wave packet 

and to a closed orbit for the lower wave packet. 



V. DISCUSSION AND CONCLUSIONS 



In this paper, we have shown that whistler waves can be efficiently trapped and advected across the magnetic field 
lines by oblique magnetosonic slow solitons with typical scale length of the order or greater than the ion inertial length. 
Oblique slow solitons carry a quasi perpendicular density perturbation that is anti-correlated to the magnetic field 
perturbation and since the propagation velocity of these solitons is much smaller than the phase velocity of whistlers, 
they can be viewed by whistlers as quasi stationary inhomogeneities. The soliton plasma density and magnetic field 
inhomogeneities then act as a true wave guide during whistler propagation. As a result, whistlers can be confined 
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in correspondence to magnetic field depletions associated to density humps (magnetic holes), as we have shown for 
frequencies uo < a;ce/2. Due to the presence of the magnetic field inhomogeneity, less strict conditions are required for 
trapping with respect to a channel provided by only a density variation. The possibility to advect whistlers marks an 
important difference with respect to the case where whistlers have been associated to non propagating mirror modes. 
Numerical results based on a two-fluid model are in good agreement with theoretical expectations. 

Slow solitons, acting as non-linear waves carrier for whistlers, provide an efficient mechanism to confine whistler 
energy in space, thus avoiding the spreading of the wave packets, and to transport the whistler energy across magnetic 
field lines at the soliton typical speed. The model we propose is related to situations often encountered in space 
plasmas and could explain multi-point Cluster spacecraft observations in the Earth's magnetosphere. The mechanism 
of whistler trapping that we have discussed relies on the "inhomogeneous, slow nature" of the wave carrier, that is, 
the plasma density and magnetic field strength inhomogeneities are anti-correlated and quasi perpendicular to the 
background magnetic field and the velocity of propagation is smaller than the whistler phase velocity. We finally 
note that other solitonic structures, propagating slowly with respect to the whistler wave packet, could in principle 
play the same role in trapping and advecting the whistlers, thus explaining space observations often showing whistler 
waves in correspondence to local minimum or maximum of the large scale magnetic field and density, respectively. 



Appendix A: Initial conditions for sIoav magnetosonic solitons 

Below the initial conditions given in the numerical code corresponding to a slow solitary wave [21^. Quantities are 
normalized to asymptotic equilibrium values outside the soliton. 



risoi 



cosh 



A/a 



VA/{12fi)x 



n=l ^Usoh Pe,i = ^Oe,i (1 ^TUsol) 



(Al) 



5^ = sin (/?o, By = cos (^o 



COS(/9o 



nsoi 



(A2) 



B. 



1 {d-^-l) Vpoi^pQ-^i) si^V^o 
Vpo de'^ (v^o - sin^ Lpo) cos ifQ 



2 V']47(12/i) tanh ^JA/(l2|l)x 



risoi 



(A3) 



Ux = VpQUsoU Uy 



VpO 



COSifo 



(A4) 



C/z ~ «i,z = -2 \/^/(12m) -^ 



(vlo-cDi^lodl-simpo^) 



cos((5o(Upo-siny)o^) 
where c^ = r(PeO + Pio) is the sound speed in normalized units and 



tanh 



^A/{12ii)x 



nsoi 



(A5) 



fJ' 



VpO (VpO - C^s) 



4de'[vlo-il+C^s)m 



{dZ^ - If sin^ ^0 
^^^ (^po - sin^ ifio) 



(A6) 



a = 



3 (fpo - cl sin^ ifo) + (f^o - sin^ y)o)[c^ + V^ (PeO + Pjo) 



(A7) 



"pO 



(l + c,^)-V(l + c2)'-4c2sinVo 



(A8) 
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Appendix B: Two-fluid cold dispersion relation 

From our two- fluid model, in the cold limit uo/k ^ Vth,e {vth,e is the electron thermal speed), we get the following 
dispersion relation for whistler waves propagating in a homogeneous magnetized plasma at rest, at an angle 6 with 
respect to the equilibrium magnetic field: 

2 (I dpm2 1 V ^^ ^ dik^?{i - cos2 ey + 2(1 + dik^)^ cos^ ^(1 + cos^ 0) + cos4 e^ \ . (bi) 



In equation (Bl) Va = B/y/n is the Alfven velocity in normalized units. 
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